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ABSTRACT: We analyze the parametric space of the constrained minimal supersymmetric stan¬ 
dard model (CMSSM) with p > 0 supplemented by a generalized asymptotic Yukawa coupling 
quasi-unification condition which yields acceptable masses for the fermions of the third family. We 
impose constraints from the cold dark matter abundance in the universe and its direct detection 
experiments, the B-physics, as well as the masses of the sparticles and the lightest neutral CP-even 
Higgs boson, m/,. We identify two distinct allowed regions with Mm > nio and mo A> Min classified 
in the hyperbolic branch of the radiative electroweak symmetry breaking. In the first region we 
obtain, approximately, 44 < tan/3 < 52, —3 < Aq/Mi /2 < 0.1, 122 < hi/,/ GeV < 127, and mass 
of the lightest sparticle in the range (0.75 — 1.43) TeV. Such heavy lightest sparticle masses can 
become consistent with the cold dark matter requirement on the lightest sparticle relic density 
thanks to neutralino-stau coannihilations. In the latter region, fixing m/, to its central value from 
the LHC, we find a wider allowed parameter space with milder electroweak-symmetry-breaking 
fine-tuning, 40 < tan/3 < 50, —11 < Ao/Mm 15 and mass of the lightest sparticle in the range 
(0.09 —1.1) TeV. This sparticle is possibly detectable by the present cold dark matter direct search 
experiments. 
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1. Introduction 


One of the most economical and predictive versions of the minimal supersymmetric standard 
model (MSSM) is the well-known - see e.g Ref. [1] - constrained MSSM (CMSSM) which employs 
just four and one half free parameters: 


signp, tan/3, M { / 2 , mo, and Aq. 


( 1 . 1 ) 


where signet is the sign of p, the mass parameter mixing the electroweak Higgs superfields Hi and 
H\ of the MSSM which couple to the up- and down-type quarks respectively, tan/3 is the ratio of the 
vacuum expectation values of Hi and H \, while the remaining symbols denote the common gaugino 
mass, scalar mass, and trilinear scalar coupling, respectively, defined at the grand unified theory (GUT) 
scale Mquj, which is determined by the unification of the gauge coupling constants. 

CMSSM can be further restricted by being embedded in a supersymmetric (SUSY) GUT with a 
gauge group containing SU (4) c and SU(2)r. This can lead to ‘asymptotic’ Yukawa unification (YU) [2], 
i.e. the exact unification of the third generation Yukawa coupling constants at the SUSY GUT scale 
Mgut- Given the experimental values of the top-quark and tau-lepton masses, the CMSSM supple¬ 
mented by the assumption of YU (which naturally restricts tan/3 ~ 50) yields unacceptable values [3] 
of the fi-quark mass for both signs of the MSSM parameter p. In Ref. [4] - see also Refs. [5-7] -, this 
problem is addressed in the context of the Pati-Salam (PS) GUT model, without the need of invoking 
departure from the CMSSM universality. We prefer to sacrifice the exact YU in favor of the universal¬ 
ity hypothesis, since we consider the latter as more economical, predictive, and easily accommodated 
within conventional SUSY GUT models. In particular, the Higgs sector of the simplest PS model [8] is 
extended so that Hi and H\ are not exclusively contained in a SU (4) c singlet, SU (2)l x SU (2)r bidou¬ 
blet superfield, but receive subdominant contributions from another bidoublet too which belongs to the 
15 representation of SU (4) c . As a result, YU is naturally violated and replaced by a set of asymptotic 
Yukawa quasi-unification conditions (YQUCs): 


/i,(Mgut) : hbiMcm) ■ /u(44gut) 


l-pa 2 /V3 


1 — p«i/\/3 


1 + \/3pcq 

y/l + \a 2 \ 2 


\ 4- |oq p 


\ + |cq p 


( 1 . 2 ) 


These conditions depend on two complex parameters cq, (Xi and one real and positive parameter p. The 
parameters cq and cxi describe the mixing of the components of the SU (4) c singlet and 15-plet Higgs 
bidoublets, while p is the ratio of then - respective Yukawa coupling constants to the fermions of the 
third family. Note that monoparametric versions of the YQUCs considered within CMSSM arising by 
taking cq = — cq for p > 0 [10] or cq = cq for p < 0 [11] are by now experimentally excluded [9,11]. 
In this talk, based on Ref. [12,13], we show that the YQUCs in Eq. ( |l.2| ) can become compatible with 
two disconnected regions of the CMSSM parameter space belonging to the hyperbolic branch (HB) of 
the radiative electroweak symmetry breaking (EWSB): The M x i 2 > mo region where the neutralino, is 
a pure bino and may (co)annihilate strongly enough with the lighter stau, f i [14,15]; and the mo A> M\ i 
region, where x acquires a sizable higgsino fraction [16-20] which enhances the x — x annihilation 
and triggers neutralino-chargino (x lx ,2 — Xfi) coannihilations. 

We begin by describing the cosmo- & phenomeno- logical requirements which we consider in our 
investigation in Sec. j2| Next (Sec. j|), we exhibit the resulting restrictions on the parameters of the 
CMSSM. Moreover, we check in Sec. ^jthe consistency with Eq. ( |L2j ) and discuss the naturalness of 
the model in Sec. [|. We summarize our conclusions in Sec. ¥\. 
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2. Phenomenological and Cosmological Constraints 


In our investigation, we integrate the two-loop renormalization group equations for the gauge 
and Yukawa coupling constants and the one-loop ones for the soft SUSY breaking parameters be¬ 
tween Mqut and a common SUSY threshold Msusy — (t\ ? are the stop mass eigenstates) 

determined consistently with the SUSY spectrum. At Msusy 5 we impose radiative EWSB, evaluate 
the SUSY spectrum employing the publicly available calculator SOFT SUSY [21], and incorporate the 
SUSY corrections to the b and r mass - see below. From M S usy to M 7 , the running of gauge and 
Yukawa coupling constants is continued using the standard model (SM) renormalization group equa¬ 
tions. The SUSY spectrum is put into micrOMEGAs [22], a publicly available code which calculates 
a number of phenomenological - see Sec. 2 A - and cosmological - see Sec. T2 - observables which 
assist us to restrict the parametric space of our model. 


2.1 Phenomenological Requirements 

2.1.1 SM Fermion Masses. After incorporating the sizable (about 20%) and less important (almost 
4%) corrections [23] to the b-quark and r-lcpton masses, we compare the masses of top-quark, m t , b- 
quark, m b and r-lepton, m T with their experimental values [24,25] 

m t (m t ) = 164.83 GeV, m b (m h ) m = 4.18 GeV, m z {M z ) = 1.748 GeV. (2.1) 

The second value is evolved up to M z using the central value a s (M z ) = 0.1185 [24] of the strong 
fine-structure constant at M z and then converted to the DR scheme with result m b (M z ) = 2.83 GeV. 

2.1.2 Collider Bounds. For our analysis, the relevant collider bounds constrain: 

• The mass m b of the lightest Higgs boson, h. The experiments ATLAS [26] and CMS [27] in the 
LHC discovered simultaneously a boson that looks very much like the expected SM Higgs boson. 
The allowed 95% confidence level (c.l.) range of m b can be estimated including a theoretical 
uncertainty of about ±1.5 GeV. This gives 

122 < m/j/GeV < 128.5. (2.2) 

• The masses of the lightest chargino, m%k, [28] and gluino, [29]: 

(a) m%± > 103.5 GeV and (b) nig >1.3 TeV. (2.3) 

2.1.3 B-Physics Constraints. SUSY contributions to observables related to fi-meson physics 
yield restrictions to the SUSY parameters. In particular, we impose the following bounds on: 

• The branching ratio BR(B S -A p 1 p ) of the process B s —>• p p [30,31] 

BR(Bj -A p + p~) < 4.2 x 10^ 9 . (2.4) 


• The branching ratio BR (b —>• 57 ) of b sy [32,33]: 

2.79 x 10 ~ 4 < BR(f? -A 57 ) < 4.07 x 10^ 4 . (2.5) 


• The ratio R (B u —> tv) of the CMSSM to the SM branching ratio of the process B u —>• tv [32] 

0.52 <R(B U ^ tv) <2.04. (2.6) 
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2.1.4 Muon Anomalous Magnetic Moment. There is a 2.9- a [34-36] discrepancy 

Sa M = (24.9 ±8.7) x 1(T 10 => 7.5 x 10 _1 ° < < 42.3 x 10~ 10 at 95% c.l. (2.7) 


between the measured value of the muon anomalous magnetic moment from its SM prediction. This 
8a^ can be attributed to SUSY contributions which have the sign of p and its absolute value decreases 
as mLsp increases. Therefore, Eq. ( |2.7[ ) hints that the sign of p has to be positive. Moreover, a lower 
[upper] bound on wlsp can be derived from the upper [lower] bound in Eq. (2/7). As it turns out, 
only the upper bound on »?lsp is relevant here. Taking into account the aforementioned computational 
instabilities and the fact that a discrepancy at the level of about 3 — a cannot firmly establish a real 
deviation from the SM value, we restrict ourselves to just mentioning at which level Eq. Q2.7[ ) is satisfied 
in the parameter space allowed by all the other constraints. 


2.2 Cold Dark Matter Considerations 


2.2.1 CDM Abundance. In the context of the CMSSM, % can be the lightest SUSY particle (LSP) 
and, thus, naturally arises as a Cold Dark Matter (CDM) candidate as long as its relic abundance does 
not exceed the upper bound on the CDM abundance deduced from the Planck satellite [37] 

£Wr< 0.125. (2.8) 


Two important mechanisms which assist to achieve Dpsp/j 2 consistent with the limit above within 
CMSSM are (i) the coannihilation of % with a particle P when a proximity between the mass /wlsp of % 
and the mass mp of P is established; (ii) the P'-pole effect which enhances the % — % pair annihilation 
procedure by an P'-pole exchange in the ^-channel when the mass of P' satisfies the relation mpi ~ 
2«?lsp- The strength of the coannihilation and the P'-pole effect processes is controlled by the relative 
mass splittings. The relevant for our cases mass splittings are defined as follows 


A p 


(j mp - m lsp) 

«?LSP 


for P = U-Xi'di and Ap> 


m P i - 2m LS p 

2w?lsp 


for P' = H, 


(2.9) 


where f\ is the lightest stau, the next-to-lightest neutralino, the lightest chargino and H the 
heavier CP-even neutral Higgs boson. The resulting (2 LS p h 2 normally decreases with these Ap’s. 


2.2.2 CDM Direct Detection. Employing the relevant routine of the micrOMEGAs package [38] 
we calculate the spin-independent (SI) and spin-dependent (SD) lightest neutralino-proton (% — p) 
scattering cross sections erf! and of®, respectively. The relevant scalar, /(-', and axial-vector, Af, form 

Ar Ar Q 

factors for light quarks in the proton (with q = u.d.s), arc taken as follows [39,40]: 

fP u = 0.018, y* = 0.026, and /£ = 0.022; (2.10a) 

A p = +0.842, \ p d = -0.427, and Af = -0.085. (2.10b) 

Data on oj* coming from large underground Xenon (LUX) experiment [41] provide strict bounds on 
the values of the free parameters of SUSY models with % owning a sizable higgsino component, as in 
our case with mo +> Mm. These data [42], however, are directly applicable in the case where the CDM 
consists of just j+s. If the jf’s constitute only a part of the CDM in the universe, the LUX experiment 
bound on the number of the scattering events is translated into a bound on the “rescaled” SI % — p elastic 
cross section where q = D^/'0.12 with 0.12 being the central value of the CDM abundance [37]. 
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3. Restrictions on the SUSY Parameters 


Imposing the requirements described in Sec. |2[ we can restrict the parameters of our model. Fol¬ 
lowing our approach in Refs. [12,13], we consider as free parameters the ones in Eq. ( |1.1[ ). The ratios 
h t /h T and hi,/h T are then fixed by using the data of Eq. ( |2. 1 [ ). These ratios must satisfy the YQUCs 
in Eq. (0) for natural values of the parameters O'], a 2 , and p - see Sec. [4j. To assure this, we restrict 
ourselves to ratios h m /h n (, m,n = t.b. t) close to unity which favor the range tan /3 > 40. We also con¬ 
centrate on the p > 0 case, given that p < 0 worsens the violation of Eq. (2.7), and scan the region 
—30 < A 0 /M 1/2 < 30. 

We localize below two separated clusters of allowed parameters categorized in the HB of the 
radiative EWSB, as justified in Sec. |3.1| : (i) The M)/ 2 > mo area at high tan/3 (43. 8 < tan/3 < 52) 


discussed in Sec. T2 and (ii) the mo Mj/ 2 area for any tan/3 in the range 40 — 50 studied in Sec. T3 


We finally - see Sec. 3A - exhibit a direct comparison of the solutions obtained in the aforementioned 
areas focusing on the characteristic value tan /3 = 48, which balances well enough between maintaining 
natural values for the h m /h n ’s and satisfying the various requirements of Sec. [5| Note that the numerical 
calculations for the soft SUSY masses become quite unstable for mo ^ M] / 2 and tan/3 > 50. 

The constraints which play an important role in delineating both allowed parameter spaces of our 
model are the lower bound on m/, in Eq. ( |2.2| ) and the CDM bound in Eq. (2.8). The M)/ 2 > m 0 region 
is further restricted by the bound on BR(B S —>• p + p _ ) in Eq. ( p.4| ) whereas the mo M) / 2 area is 
additionally bounded by the LUX data and the limits on m%± and nig in Eq. ([U3|a) and (b). In the 
parameter space allowed by these requirements, all the other restrictions of Sec. [| are automatically 
satisfied - with the exception of the lower bound on 8a u in Eq. (2.7). 


3.1 Elliptic Versus Hyperbolic Branch 

The classification - see e.g. [19,20] - of the various solutions of the radiative EWSB condition is 
based on the expansion of p 2 in terms of the soft SUSY breaking parameters of the CMSSM included 
in Eq. ([0]). Indeed, using fitting techniques, we can verify the following formula 

p“ +M ^/2 — cothq + Ci/ 2 M 2 ^ 2 + c am^()3Tj/ 2 , (3.1) 

where the coefficients co,Ci/ 2 ,ca, and cam depend basically on tan/3 and the masses of the fermions of 
the third generation. These coefficients are computed at the scale Msusy - see Sec. ^ and, therefore, 
inherit a mild dependence on the SUSY spectrum too. From Eq. ( |3.1p , we can easily infer that the 
SUSY breaking parameters are bounded above for fixed p, when the quadratic form in the right-hand 
side of this equation is positive definite. This is the so-called ellipsoidal branch (EB) which is highly 
depleted [20] after the discovery of li with m/, in the range of Eq. (^). On the other hand, in the HB 
region favored by Eq. ( |2.2[ ), co is negative and, consequently, mo can become very large together with a 
combination of Ao and M { / 2 with all the other parameters being fixed. In this case the soft parameters 
lie on focal curves or surfaces. Near the boundary between the EB and HB regions, co is very close 
to zero and, thus, mo can become very large with all the other parameters fixed. These are the so- 
called focal points. Moreover, there is a region where the soft SUSY breaking mass-squared m 2 Hi of Hi 
becomes independent of the asymptotic value of the parameter mo- This is called the focus point (FP) 
region [18]. In the large tan/3 regime under consideration, we have [1] mjp ~ —p 2 — Af|/2 and so no 
distinction between focal and focus points can be established. 
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Region 

co 

c l/2 

ca 

cam 

Co/10 8 

C A / 10 7 

Mi/2 > m Q 

-0.059 

0.955 

0.102 

-0.277 

-1.05 

6.14 

mo 2> Mi/ 2 

-0.0689 

0.848 

0.0982 

-0.254 

33.3 

-2.34 


Table 1: The c’s in Eq. (3.1) and C’s in Eq. (3.2) for tan /3 = 48 and M , /2 > otq or mo M l / 2 . 


To get an idea of how our solutions presented in Secs. [T^-|3~4| are classified into these categories, 
we display in Table [I] the values of the coefficients c in Eq. ( |3. 1 1 ) for the two representative cases of 
Table ^ corresponding to tan/3 = 48 - see Sec. 3.4. We obtain Co < 0 in both cases and so we expect 
that our solutions belong to the HB region. To illustrate the emergence of the relevant focal curves we 
diagonalize the quadratic form in the right-hand side of Eq. ( |3. ip keeping, e.g., M x / 2 fixed and using 
mu and Ao as varying parameters. Then, Eq. ([371]) can be cast in the following form 


mllC 0 +A 2 /C a = I 


(3.2) 


where Aq = A + camMi/ 2 /2ca, Cq = JU 2 /co, and Ca = JU 2 /ca with ft 2 = ft 2 +M|/2 — C]/ 2 M 2 p and 
c | /2 = C\j 2 — c 2 am /Aca- The numerical values of Co and Ca also listed in Table [I]. Since in both cases, 


one from Co and Ca is negative, mo and Ao can vary along an hyperbola parameterized by Eq. ( |3.2| ) and 
so both (M | / 2 > mo and mo A> M\ / 2 ) our solutions in Table |] belong to the HB region. 

3.2 Mi / 2 > m 0 Region 

Initially, we concentrate on the M\ /2 > n?o region and delineate in the left plot of Fig. [1] the allowed 
(shaded) areas in the M l / 2 — mo plane for tan/3 = 48 and various Ao/M l / 2 ’s indicated therein. The 
lower boundaries of these areas corresponds to Af, = 0; the areas below these boundaries are excluded 
because the LSP is the charged tj. The upper boundaries of the areas come from the CDM bound in 
Eq. (243), while the left one originates from the limit on BR(/? V —>• /./ 1 /i~) in Eq. ([2.4|). The upper 


right corners of the areas coincide with the intersections of the lines = 0 and £2 LSP /r = 0.125. We 
observe that the allowed area, starting from being just a point at a value of Ao/M l / 2 slightly bigger than 
—0.9, gradually expands as Aq/Mi/ 2 decreases and reaches its maximal size around Ao/M| 2 = —1.6. 
For smaller Ao/Mj^’s, it shrinks very quickly and disappears just after Ao/M)/ 2 = —1.62. The fact that 
the allowed regions are narrow strips along the lines with A 2| = 0 indicates that the main mechanism 
which reduces £2 lsp It 2 below 0.125 is the % — i\ coannihilations. Namely, the dominant processes 
are the f 2 f 2 coannihilations to bb and tt which are enhanced by the ^-channel exchange of H. with 
Ah ~ 1.1 - see also Table [2. 

Extending our analysis to various tan /3 ’s for A 2| = 0 we can obtain a more spherical view of the 
overall allowed region of the model for M\ / 2 > mo- This is because A^ = 0 ensures the maximal 
possible reduction of Tlisph 2 due to the — f| coannihilation and so we find the maximal M\ / 2 or 


/jjlsp allowed by Eq. (243) for a given value of Ao/M\/ 2 . The relevant allowed (hatched) regions in the 


Mi /2 — Ao /Mi /2 plane are displayed in the right plot of Fig. [l[ The right boundaries of the allowed 
regions correspond to £2 lsp It 2 = 0.125, while the left ones saturate the bound on BR (R s ) in 

Eq. ( |2.4[ ). The almost horizontal upper boundaries correspond to the sudden shrinking of the allowed 
areas which is due to the weakening of the //-pole effect as Ao/M } i 2 drops below a certain value 
for each tan/3. The lower left boundary of the areas for tan/3 = 44, 45, and 46 comes for the lower 
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M 1/2 (TeV) 



FIGURE 1: Allowed (shaded) areas in the Min — mo [Min — Aq/Mi n] plane (left [right] plot) for Min > mo- 
In the left [right] plot we use tan/3 = 48 [Af, = 0] and various Aq/Mi n’s /tan/3 ’s] indicated in the graph. The 
dot-dashed [dashed] line corresponds to mh = 125 [126] GeV. 


bound on m/, in Eq. (2.2), while the somewhat curved, almost horizontal, part of the lower boundary 


of the area for tan/3 = 44 originates from Eq. (2.8). The dot-dashed and dashed lines correspond to 
m/j = 125 and 126 GeV respectively. We see that the nth’s which are favored by LHC can be readily 
obtained for 47 < tan/3 < 50. In the overall allowed region we obtain 122 < m/j/GeV < 127.23 and 
746.5 < ntLsp/GeV < 1433. Also ~ (0.35 — 2.76) x 10~ 1() and so, Eq. ( |2.7[ ) is satisfied only at the 
level of 2.55 to 2.82 — a. 


3.3 m 0 >Mj/ 2 Region 

The interplay of the various requirements of Sec. ^]in the mo 2> M) ; / 2 region can be easily under¬ 
stood from Fig. [2| where we present the (shaded) strips in the Mm — mo plane allowed by Eqs. 


- (±8) for tan/3 = 48 and several Aq/M\/ 2 s indicated in the graph - note that no restrictions from 
LUX data are applied to this plot. The upper [lower] boundary along each of these allowed strips arises 


from the limit on m%± [Glsp/* 2 ] in Eq. 


|a) [Eq. ( |2.8| )1. On the other hand, the lower limit on m/, in 
Eq. ( |2.2[ ) causes the termination of the strips at low values of m 0 and Mp 2 . whereas their termination 
at high values of mo is put in by hand in order to avoid shifting the SUSY masses to very large values. 
The solid lines indicate solutions with m/, = 125.5 GeV - see Eq. <[2.2|). From this figure, we easily see 
the main features of the mo A> M { / 2 region: mo spans a huge range (4 — 15) TeV, whereas M { / 2 \[l ] 
remains relatively low (1 — 6) TeV [(0.1 — 1) TeV]. We observe also that as Aq/Mi/ 2 increases from 
—2 to 2 the allowed strip moves to larger M|, 2 ’s and becomes less steep. 

Varying continuously Aq/M\/ 2 for tan/3 = 48, m/, = 125.5 GeV and taking into account the LUX 
data, we depict in the right plot of Fig. j2j the overall allowed region of the model for mo A> Mm 
in the mLSP — Aq/Mi/ 2 plane. On the solid and dashed line, the bounds on m%+ in Eq. (]2.3|a) and 
on G[ si'/; 2 in Eq. ( p.8| ) are saturated, whereas the restriction from the LUX data on q yields the 
dotted boundary line. Finally, the double-dot dashed boundary lines from the limit on trig in Eq. (|T3|b) 
provide the maximal and minimal Ao/M^’s. Note that the allowed regions are almost symmetric about 
Aq/Mi/ 2 ~ 2.5. Also, we remark that // remains almost constant ~ 100 ±20 GeV on the solid lines 
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FIGURE 2: Allowed (shaded) areas in the M l n — mo /wlsp — Aq/Mx/i] plane (left [right] plot) for tan/3 = 48 
and mo M\ / 2 . In the left plot we take various Aq/Mi/j’s indicated in the graph and the black lines correspond 
to nih = 125.5 GeV. In the right plot we set mi, = 125.5 GeV. 


from Eq. ( J2.3[ ), while it reaches about 1 TeV when the bound in Eq. ( |2.8[ ) is saturated. Close to the 
latter portion of the parameter space £ 2 lsp h 2 calculation is dominated by the x/X 2 ~ X\ coannihilation 
processes whereas in the region where m% is near its lower limit in Eq. <[2.3|a), the % — % annihilation 
processes contribute more efficiently to the resulting kl^sph 2 . All in all, we obtain —12. 8 < Ao/My/2 < 
15.8 and 92 < n?LSp/GeV < 1084.2. In this area 8 a^ is well below the lower limit in Eq. (p~7|), i.e. 


8an ~ (0.04 — 0.27) x 10 . Therefore, Eq. (2.7) is satisfied only at the level of 2.83 to 2.86 — a. 

si 
XP 


In the mo A> M x / 2 region, where % has a significant higgsino component - see Table § -, a: 


is 


dominated by the /-channel Higgs-boson-exchange diagram contributing to the neutralino-quark elastic 
scattering process - for the relevant tree-level interaction terms see e.g. the appendix of Ref. [ 6 ]. Espe¬ 
cially for large tan/3’s, which is the case here, the couplings of H to down-type quarks are proportional 
to tan/3 and so are the dominant ones. More explicitly, ct|* behaves as 


j§i , 
J XP 


tan 2 p\N U i\ 2 \N h3 \ 2 /m 4 H , 


(3.3) 


where Aqi, N\o, and ( 3 4 ) arc the elements of the matrix N which diagonalizes the neutralino mass 
matrix and express the bino, wino, and higgsino component of %, respectively. As a consequence, ct|* 
can be rather enhanced compared to its value in the > «?o region, where % is a pure bino. 

This conclusion can be clearly induced by Fig. |3| where we show the allowed (shaded) regions in 
the »?lsp — q cr|p plane for tan/3 = 48. The left [right] panel corresponds to Aq/A/j / 2 < 0 [Ao/M 1/2 > 0]. 
The numbers on the various points of each boundary line indicate the corresponding values of Aq/Mi/ 2 . 
The (black) solid, double dot-dashed, and dashed lines correspond to the bounds from Eqs. (p~3|a), 
d2.3| j), and ( |2.8p . respectively - cf. right plot in Fig. The dotted lines arise from the LUX data, 
whereas the dot-dashed lines give the lowest possible q (jA in each case, with q ranging from from 
about 0.013 to 1 (along the dashed line). From these graphs, we infer that, the maximal [minimal] 
is located in the upper right [lowest left] corner of the allowed regions, at the junction point of 
the dashed and dotted [dot-dashed and solid] lines. The overall minimum of 8, o-, 1 is obtained in the 
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—<— nu = 103.5 GeV 

x- 

Q LS P h2 = 0-125 

- rru= 1.3 TeV 

■ minimal ^a sl 

- - - LUX Data 

XP 


FIGURE 3: Allowed (shaded) regions in them^sp ^< 7 ?^ plane for mo Mi/ 2 , m/, = 125.5 GeV and tan /3 = 48. 
The left [right] panel corresponds to Aq/M\ / 2 < 0 [Ao/Min > OJ and the values ofAo/M \/2 at the various points 
of the boundary lines are indicated. The conventions adopted for the various lines are also shown. 

right plot of Fig. |3] whereas the maximal one is practically the same in both plots. Namely we obtain 

1.66 x 1(T 12 < ^crfp/pb < 1.03 x 1(T 8 . (3.4) 

The obtained values of <Ejo|* are within the reach of the forthcoming experiments like XENON1T [43] 
and LUX-ZEPLIN [44], whose the planned sensitivities are depicted in Fig. |j by a gray solid and 
dashed line respectively. Note, finally, that similar results, as regards the allowed ranges of »?lsp and 
Ao /Mi /2 and the predicted q a^ l p can be obtained [13] for other tan /3's in the range 40 — 50 too. 

3.4 Mi / 2 > mo Versus m 0 > Af 1//2 Region 

Comparing the right panels of Figs. [I] and || we notice that the allowed areas in the M\ / 2 > m o and 
mo 2> Mi /2 regions share common Ao/M^’s for tan/3 = 48. Focusing on a such Aq/Mi / 2 value, e.g. 
Aq/Mi /2 = —1.5, and selecting the remaining input parameters so that we achieve the central values 
Glsp^ 2 = 012 and m/, = 125.5 GeV, we can implement a direct comparison between the M ] / 2 > mo 
and mo 2> M x / 2 solutions. This is done in Table [3. where we arrange the values of the input and some 
output parameters, the mass spectra and some low energy observables for two characteristic points of 
the allowed parameter space our model. The various masses of the SUSY particles (gauginos/higgsinos 
X, X 2 , X 3 , X 4 , % 2 ’ 8’ squarks t u h, b\, b 2 , u L , Hr, d L , d R , and sleptons fi, f 2 , V T , e L , e R ) and the 

Higgs particles ( h , H, H ± , A) are given in TeV - note that we consider the first two generations of 
squarks and sleptons as degenerate. From the values of the various observable quantities, we can verify 
that all the relevant constraints, but the one of Eq. ([T7|), are met - cf. Sec. |[ For the interpretation of 
our results, mainly on £2 lsp^ 2 > w e also list the values of the various A/2s in Eq. ( ]2.9[ ), the bino, (Ah 1 1 2 , 
and the higgsino, Vi a) 2 + |/V| 4 | 2 , purity of %. We also include an estimate for the EWSB fine-tuning 
parameter Aew - see Sec. |5| 
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Input Parameters 

tan/3 

48 


—Aq/Mi/2 

1.5 


M 1/2 /TeV 

2.821 

2.157 

m o /TeV 

2.522 

9.219 

Output Parameters 

ht/h t(Mgut) 

1.117 


h-b / hr (A/gut ) 

0.623 


h t /hb{M gut ) 

1.792 

1.45 

Id /TeV 

3.514 

0.936 

Aew 

2972 

216.2 

(%) 


615 

(%) 

‘ 1 - KafFH 

94 

V (») 

89.96 

1.56 

Afc (%) 

89.95 

1.928 

|M,i| 2 

1 

0.145 

|M, 3 | 2 + |M,4| 2 

0 

0.852 

Masses in TeV of Sparticles and Higgses 

<v / v ^ 

A 5 A 2 

1.305,2.478 

0.943,0.960 

y0 yO 

A3 5 A4 

3.509,3.512 

1.034,1.959 


2.479,3.512 

0.956, 1.959 

g 

6 

4.936 

h,h 

3.998,4.729 

6.309,7.270 

bi,b 2 

4.692,4.772 

7.267,7.887 

ml, «r 

5.880,5.625 

10.1,10.004 

dL,dR 

5.880,5.592 

10.1,9.992 

Ti,t 2 

1.309,2.661 

6.749,8.202 

CL,eR 

3.162,2.744 

9.359,9.276 

V T , V e 

2.656,3.160 

8.201,9.359 

h,H 

0.1255,2.640 

0.1255,3.67 


2.641,2.640 

3.671,3.67 

Low Energy Observables 

10 4 BR(Z> —»• sy ) 

3.27 

3.3 

10 9 BR (B s ->ii + h~) 

3.74 

3.02 

R (B u -A tv) 

0.984 

0.991 

io 10 ^^ 

0.68 

0.227 

Olsp h 2 

0.12 

°f P / P b 

3.35 • 10 -12 

7.28 • 10 -9 

a?D/p b 

6.67- IQ’ 10 

8.57 • 10~ 6 


TABLE 2: Input/output parameters, sparticle and Higgs masses, and low energy observables for tan j3 = 48, 
Aq/Mi /2 = — 1.5 and M x / 2 > hiq (second column) or mo M\/ 2 (third column). 
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The results of the two columns of Table |2] reveal that the M \ ; / 2 > mo and mo 3> M\ /2 solutions ex¬ 
hibit a number of important differences: First of all, in the mo 3> Mj/ 2 region, mo acquires considerably 
larger values, while ft remains quite smaller than its value for M l / 2 > mo- As a consequence, the Fliggs 
bosons (//, H ± , and A) acquire larger masses and the whole sparticle spectrum, with the exception 
of the neutralinos and charginos, becomes heavier. As a by-product, the various observables besides 
ffi sr>/i 2 acquire values closer to non SUSY ones. E.g., tkq, is even smaller than in the M\/ 2 > mo 
region. Note that the latter region is tightly constrained by BR(B S —» ft ft ), which is well suppressed 
for mo 2 > M l / 1 . Similar values for the ratios h m /h n with m.n = t.b.x are obtained in both cases with 
results slightly closer to unity for mo M { / 2 . 

To shed more light on the mechanisms which ensure Q.^sph 2 com¬ 
patible with Eq. ( |Oi| ) for M\ / 2 > mo and mo S> M\ / 2 , we amange in Ta¬ 
ble |] the relative contributions beyond 5% to 1 /TUsp/; 2 of the various 
(co)annihilation processes for the inputs of Table |2| From these results, 
we infer that a synergy between % ~ T 1 coannihilation and the //-funnel 
mechanism is well established in the M\ / 2 > mo region thanks to the quite 
suppressed Aq and A h achieved. On the other hand, for mo M\p_ the 
X/X 2 — X\ coannihilations are activated because of the low A^+’s and 
A^’s. However, in the latter case the annihilation channels to IT IT h 
and ZZ conserve their importance due to the large higgsino mixing of %. 
Due to this fact /7 ?lsp is confined close to /./ and, as a bonus, the result¬ 
ing cr|p’s are accessible to the forthcoming experiments [43,44]. On the 
other hand, for A/) /2 > mo, X is an almost pure bino with mLsp — 4/] / 2 /2 
and rr|* well below the sensitivity of any planned experiment. In both 
regions, finally, turns out to be much lower than the reach of Ice- 
Cube [45] (assuming x ~ X annihilation into IT IT U and the expected 
limit from the large DMTPC detector [42]. Therefore, the LSPs predicted 
by our model can be detectable only by the future experiments which will 
release data on <7^ l p and for mo 3> Mi/ 2 . 

4. Deviation from Yukawa Unification 

In the overall allowed parameter space of our model [12,13], we find the following ranges for the 
ratios h m /h n with m.n = l.b. t: 

0-98 < h t /h T < 1.29, 0.60 < h b /h T < 0.65, and 1.62 < h t /li b < 2 for M l/2 > m 0 ; (4.1a) 
1.00 < li t jhx < 1.50, 0.75 < h b /h x < 0.79, and 1.20 < h t /h b < 2 for m 0 >M 1/2 (4.1b) 

We observe that, the required deviation from YU is not so small and turns out to be comparable to the 
one obtained in the monoparametric case - cf. Ref. [10]. In spite of this, the restrictions from YU 
are not completely lost but only somewhat weakened. Actually, our model is much closer to YU than 
generic models with lower tan /3 ’s where the Yukawa coupling constants can differ even by orders of 
magnitude. Also, the deviation from YU is generated by Eq. ( |L2| ) in a natural, systematic, controlled 
and well-motivated way. 


/ 2 > m 0 Region 

Ti -ff —f bb 

69% 

fl f * —f XX 

15% 

mo >M 1/2 Region 

XX\ ud 

18% 

Xixt ->• Ud 

8% 

xxt -> tb 

7% 

XXt ->• 

6% 

XX^W-W + 

6% 

XX 

5% 


TABLE 3: Processes which 
contribute to I/Olsp^ 2 more 
than 5% and their relative 
contributions for M { / 2 > mo 
and mo ^ Mp 2 . 
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8 2.0 

E 

1.5 


1.0 



-1.5 -1.0 


-0.5 0.0 

Re a _ 


0.5 1.0 



3.0 


2.5 


8 2.0 

E 

1.5 


1.0 


-1.5 


A p = 0.5 
° p = 0.4 
D p = 0.3 


- 1.0 


-0.5 0.0 

Re a. 


0.5 


1.0 



FIGURE 4: The complex parameters cq and a 2 for various p’s indicated in the graphs for tanf) = 48, Aq/Mi /2 = —1.5, 
m/, = 125.5 GeV and Mj/ 2 > 'no (upper plots) and mo 2> Mf/ 2 (lower plots) - see Table 2. 


To show explicitly it, we below extract values of p, cq, and oc 2 , which leads to the ratios h m /h n with 
m,n = t.b. T, encountered in the characteristic examples presented in Table |2[ Since from Eq. ( |1.2| ) we 
have only two equations and five real unknowns we can find infinitely many solutions. Some of these 
solutions are shown Fig. [| for Mm > mo (upper plots) and mo M)/ 2 (lower plots) and various p’s 
indicated therein. Since the equation for lib/h x depends only on the combination pfX\ its solutions, for 
fixed p, lie on a certain curve in the cq complex plane, as shown in the left upper and lower panels of 
Fig. §. For each cq and p in these panels, we find various a 2 ’s, depicted in the right panels of Fig. f|, 
solving the equation for h t /h x . Observe that the equation for h t /li r depends separately on a 2 an d p 
and, thus, its solutions do not follow any specific pattern in the a 2 complex plane. Scanning the range 
of p from 0.3 to 3 and we can find solutions in the cq and a 2 planes only for the lower values of this 
parameter (up to about 0.6) for both allowed regions of the model. These solutions are very similar to 
the ones displayed in Fig. j| for all the possible values of the ratios of h m /h n with rn.n = t.b. T allowed 
by the constraints of Sec. |2. Consequently, we can safely conclude that these ratios can be readily 
obtained by a multitude of natural choices of the parameters p, a 1 , and a 2 everywhere in the overall 
allowed parameter space of the model. 
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1 1.1 1.2 1.3 0.85 0.9 0.95 1 1.05 1.1 


itilsp (TeV) m L sp (TeV) 

FIGURE 5: M\/ 2 and Agw as functions of /mlsp f° r m h = 125.5 GeV, £2 lsp^ 2 =0.125 and various A o /M\ y 2 ’ s indicated on 
the curves for the Mj / 2 > mo (left panel) and the mo 2> Mj / 2 (right panelj region of the model. 


5. Naturalness of the EWSB 

The fact that, in our model, Mw 2 , mo, and /./ generally turn out to be of the order of a few TeV 
puts under some stress the naturalness of the radiative EWSB giving rise to the so-called little hierarchy 
problem. To quantified somehow this issue we introduce the EWSB fine-tuning parameter 

/ |C;| \ , \ ( o m H, ml, tan 2 /3\ 

4EW ^ m "(Mp) With (5-1) 

Here i = fl.H\.H 2 and »?//, is the soft SUSY breaking mass of Hj with j = 1,2. In most of the parameter 
space explored, Aew is dominated by the term C u . 

Focusing on the values of the parameters which ensure Gpsp/; 2 ~ 0.125, we present, in the left 
[right] panel of Fig. |5| M\j 2 (solid line) and Aew (dashed line) as functions of »jlsp for Mm > n?o 
[mo A> AY| , 2 !, tan/3 = 48, m/, = 125.5 GeV, and negative Ao/M^’s indicated in the graphs. The 
Aq/M 1 / 2 , s fulfilling these conditions are shown in Fig. 3 of Ref. [12] for M l/2 > mo and in the left panel 
of Fig. |3] for mo iVf\/ 2 . We clearly see that, in both regions, the resulting Aew is almost constant and 
Aew ~ 2000 for M) ,- 2 > mo whereas Aew ~ 200 for mo A> Mj n in agreement with the values shown in 
Table |2| In other words, Aew in the M u2 > mo area becomes about a factor of ten larger than its value 
in the mo M) j 2 area despite the fact that the resulting /ulsp’s are comparable. The crucial difference 
between the two regions, though, is the lower ft’s encountered for mo ($> M { / 2 - see Table |] - which 
largely influences the Aew computation - see Eq. ( |5.1| ). We can conclude, therefore, that the mo S> A/, / 2 
solutions are more natural regarding the EWSB fine-tuning than those for M\/ 2 > mo- 

6. Conclusions 

We investigated the compatibility of the generalized asymptotic YQUCs in Eq. (|P|), which yield 
acceptable masses for the fermions of the third family, with the CMSSM for ju >0 and 40 < tan /3 < 50. 
We imposed phenomenological constraints originating from the mass of the lightest neutral CP-even 
Higgs boson, the lower bounds on the masses of the sparticles, and B-physics. We also considered 
cosmological constraints coming from Lh sp/? 2 and the LUX data on q 
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We found that % can act as a CDM candidate in the following two separated regions classified in 
the HB of the EWSB: 

• The M 1/2 > n?o region, where the LSP turns out to be an essentially pure bino and Qi^ph 2 is 
reduced efficiently, thanks to //-pole enhanced stau-antistau coannihilations, so that it is compat¬ 
ible with the recent data on BR (B s —>• /./ 1 fl~). The LHC preferred values /??/, ~ (125 — 126) GeV 
can be accommodated for 48 < tan/3 < 50, whereas »?lsp comes out to be large (~ 1 TeV). As a 
consequence, the % direct detectability is very difficult and the EWSB fine-tuning becomes rather 
aggressive since Aew 1 ~ 0.035%. 

• The mo i$> M\ ,- 2 region, where the LSP is a bino-higgsino admixture and has an acceptable 
Q.I SPIr thanks to % — % annihilations (for low «?lsp's) and the x/Xi~ %\ coannihilations (for 
large /wlsp’s). Fixing /??/, = 125.5 GeV favored by the LHC, we found a wider allowed parameter 
space with 40 < tan/3 < 50, —11 < Ao/M x /2 < 15, 0.09 < n?LSp/TeV <1.1 and milder EWSB 
fine-tuning since Aew 1 ~ 0.5%. The LSP is possibly detectable in the planned CDM direct 
search experiments which look for a|*. 

In both cases above, the restriction on 8a^ is only satisfied at a level of above 2 — a and the 
required deviation from YU can be easily attributed to a multitude of natural values of the relevant 
parameters within a PS SUSY GUT model. 
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